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' Abstract. The Taylor coefficients c and 

d of the Pion EM form factor are con- 

' strained using analyticity, knowledge of 

the phase of the form factor in the time- 
like region, 4m^ < t < tii^ and its value 
at one space-like point, using as input the 
{g — 2) of the muon. This is achieved us- 
ing the technique of Lagrange multipli- 
ers, that gives a transparent expression 

' for the corresponding bounds. We present 

a detailed study of the sensitivity of the 

' bounds to the choice of time-like phase 

and errors present in the space-like data, 

, taken from recent experiments. We find 

that our results constrain c stringently. 
We compare our results with those in lit- 
erature and find agreement with the chi- 
ral perturbation theory results for c. We 
obtain d ~ 0(10) GeV"'' when c is set to 
the chiral perturbation theory values. 

1 Introduction 

The Pion form factor continues to be of current inter- 
est [T1[51[3J|3] . In [2] , we developed a framework for obtain- 
ing constraints on the low-energy expansion coefficients of 
the Pion form factor using data from the space-like region 
{t < 0) l5,,,6i7,8jiS], suitably extending an earlier work of 
Raina and Singh [lOj . It was shown using arguments of an- 
alyticity for the Pion form factor and a reliable estimate 
of the pionic contribution to the {g — 2) of the muon, that 
it is possible to isolate allowed regions in the c — d plane, 
where c and d arc the Taylor coefficients in the low-energy 
expansion of the Pion EM form factor given as, 

F^it) = F^iO) + l{rl)t + ct^+dt^ + --- (1) 



The form factor is analytic everywhere in the complex t 
plane except for the cut along (4mJ,oo). Caprini [TTj has 
employed the phase of the time-like data along a part of 
the cut and the QCD polarization function iT' in order 
to constrain the coefficients (c, d) . Each of the above in- 
dependently use either pure space-like data or phase of 
time-like data, employing the method of Lagrange multi- 
pliers. In this work, we consider the problem of simulta- 
neous inclusion of phase of time-like data and one space- 
like datum. The constraints are introduced through La- 
grange multipliers: the constraints from the phase of time- 
like data is introduced through an Omnes function, while 
the space-like constraint has a simple linear form as seen 
in [2]. Caprini also advances a method of implementing 
data from time-like phase as well as modulus that does 
not use the technique of Lagrange multipliers in . 

As shown in pl[TT| , the coefficients constrained by the 
normalization of i<V and the value of the pion charge ra- 
dius satisfy the following inequality 

oo 

i^j2'^c,,r<i. (2) 

n=0 

The c„'s are functions of the Taylor cocfhcients, c, d, and 
an "outer function" related to the observable of interest, 
viz., {g — 2) of the muon, 77' etc. The inequality in Eq. ([2]), 
is modified by the additional constraints provided by the 
space-like data and the phase of time-like data along a 
part of the cut. We already saw in ref. [2], the modifica- 
tion to Eq. ([2]) that comes through space-like constraints 
alone, where the constraints are introduced through La- 
grange multipliers at. Eliminating them results in a simple 
determinantal equation for the bounds. Here, we also take 
into account the phase of the form factor 5{t), for t in the 
range tj, to tin, where tin = 0.8 GeV^. The Fermi- Watson 
final state theorem implies that in the elastic region the 
phase of the form factor is given by the P-wave tt — tt 
scattering phase shift Si . For a detailed discussion of the 
theorem we refer to ref. [1] . Analogous to [11] , we vary the 
coefficients for n > -I- 1, assuming that the first A^ co- 
efficients are known through the normalization of and 
the value of r^. We express the resulting constraints in a 
form where improvements, input by input, are manifest. 
Setting A" = 3 and using an input for / yields the ellipse 
in the c — d plane. 

A good check for our results is obtained by setting 
N = 2 and comparing with the treatment of Bourrely and 
Caprini |12j , who have considered a problem in the tt K 
sector and have included the tt K phase shift and con- 
straints from the Callan-Treiman point, which is a time- 
like datum in the analyticity domain. 

This paper has been divided into the following sec- 
tions. In Section [21 we describe the Lagrange multiplier 
technique and obtain expressions for the constraints and 
bounds. In Section [3l we discuss the model of the phase 
of time-like data we choose and the data sets we employ 
for the space- like part. We present our results in Section [Jj 
and summarize in Section [5l 
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2 Formalism 

The pion contribution to the niuon anomaly is given by: 



a^{n+TT-) = - I dtp{t)\F^{t)\ 



(3) 



where <^ = 4mJ is the branch point of the pion form factor 
and 

M = ^'^^m > (4, 



where, 



K{t)= I du(l-u)u'^(l-u+ —!^)-\ (5) 







t 



We use the following map from the i-plane pillOj. that is 
cut from along the real t axis, to the complex z-plane 
(region \z\ < 1), 



z — 1 / f — 
= i\ 



2 + 1 



(6) 



This map takes the point t^r to 1 and the point at infin- 
ity to —1, as seen in Fig. (jT]). Using this map and the 
definitions: 

/(^) = F^it) (7) 

p{z) = Pit) (8) 

the pionic contribution to the muon anomaly can be writ- 
ten as 



1 



2n 



S\\2 



{n+n-) = — d0w{9)\fie^')\ 



2tt Jo 
where, 

w{e) = Ami sec2(6'/2) tan(6l/2) p {e'") > 0. 
We now consider a function h{z) defined as: 
H^) = f{z)wT,{z) 

where. 



Wti{z) = exp 



1 r^"" .„e*' 



47r Jo e" 



■Itiw{6) 



Then Eq. ([9]) can be written as: 



ap(7r+7r ) = 



1 

2^ 



dd\h{e''^)\^^ 



(9) 

(10) 
(11) 

(12) 
(13) 



Now h{z) is analytic within the unit circle |z| < 1 and 
for real z, h(z) is real. Therefore h{z) can be expanded as 
follows: 

h{z) = flQ + aiz + a2Z^ + ■ ■ ■ (14) 

where ag, ai • • • are real coefficients. Therefore, in the an- 
alytic region, \z\ < 1, a^(7r"'"7r~) can be written as: 



a^(7r+7r ) = + of + ■ 



(15) 



a consequence of Parscval's theorem of Fourier analysis. 
The expansion coefficients a„ can be obtained from a Tay- 
lor expansion of the function h(z) in terms of f{z) and 
Wjr^z). The coefficients are given by 



ao = h{0) = w^(0), 
2 



(16) 



ai^h'iO) = wUO) + ^rlUw^{0), (17) 



02 = 



h"{0) 



2! 



w-(0) ( --r% + 32ctl 



and 



0-3 



h'"{0) 



2<(0) ( ^rX 



<(0) 



(18) 



3! 



[w,{0) {I2r% - 384 c4 + 384^4)] 



3<(0) i^~-r% + 32ctl 



Here, the expansion coefficients satisfy 



oo 

E 

Tt = 



[an? 



(19) 



(20) 



Given /. Eq. ([^0]) yields constraints on the expansion co- 
efficients of the form factor. Including up to the second 
(third) derivative for F.^{t) results in constraints for c (c 
and d). Now if we divide Eq. (^0]) by /, then, 



where, 



(21) 



(22) 



Any c„ which satisfies /^g < 1 is allowed and the equality 
gives the bound. We have already seen that this yields an 
ellipse in the c — d plane plfTT] . 

The time-like phase is the argument of the form factor 
along a part of the cut. Let us assume that we know the 
phase in the region t^^ <t < ti^. We assume that. 



(23) 



in accordance with Watson's theorem, where Sl{t) is the 
two-pion elastic scattering phase shift and ti^ = 0.8 GeV^ 
[TT|[T3]. The z map, as defined in Eq. ([6]), maps t^^ 1 
^9^0 and Un ^ z-^ ^ ^ 0-^, where z = e'^ The 
upper part of the cut is mapped on to the upper half circle 
and the lower part of the cut to the lower half circle. As 
a result, 

Avg [F^ire'')] = 6l{9), < 6 < 6,,, (24) 
Arg [F^ire^')] = -Sl{9), {2n - (?;„) < 9 < 2tt (25) 
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t-plane 
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time-like 



upper-lip 
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Fig. 1. z — t map for z — e'". The broken blue arc represents the time-like region, the dashed-dot red line the space-like region 
and the dashed-dot-dot-dashed green line the time-like region along positive real t axis. The solid blue arc highlights the region 
in the z plane which contributes to the phase data used here. 



as shown in Fig. ((T|). 

The phase is introduced through the Omnes function, 
denoted as Ot^{z) in the z-plane, by. 



Ojriz) ~ exp 



1 — ze' 



(26) 



where 



= -5l if 2n - 6ii„ <e <2Tr. 



(27) 



Since the phase of the form factor along the cut (i.e. t^r < 
t < tin) is compensated by the phase of the two pion 
scattering phase shifts, the following condition holds. 



Im lim 



1 



= 



(28) 



which is the constraint equation coming from the time-like 
phase of the form factor. 

In order to include information on phase of time-like 
data as well as space-like data, we consider the following 
Lagrangian 



C 



2 ^ 

1 °° 



c„ lim 

n=0 



A(6l) I M^(6l) |Im[[PT/(6l)] " V"e™''](i6l 



a{J ~ ^ c„z"), 



n=0 



where, 
and 



w{0) = w{C) = wAOOAC) 



(29) 



(30) 



C = cxp(i6i). (31) 

J = h{z)/^/l is the space-like datum mapped on to the z 
plane as defined in Eq. ([6]) and / is the bound from either 



or 77'. We note that the space-like region is mapped 
on to the negative real axis as can be seen in Fig. ([T|). 
The value of z and h{z) used is given in Table [TJ As in [21 
[TT| . we use a value of / = 75 x 10~^ in the following. As 
mentioned in [2l[T0]. the space- like constraints are linear 
and can be written as: 



(32) 



ri=0 



and time-like constraints arc given through the following 
equation 



T1=0 



iJm 



lim [[1^(61)]- 







(33) 



We minimize the Lagrangian in Eq. ([M)) for ri > -f- 1 
where N represents the terms already constrained using 
the normalization of TV and the pion charge radius r,r- 
In other words, we first constrain the coefficients though 
Eq. (PT|) . truncating it at N and then impose additional 
constraints through the time- like phase and space- like data. 
We have already seen that for A^ = 3, this yields an el- 
lipse in the c — d plane pilTT] . Denoting the coefficients for 
n > A^ -I- 1 as c„ we have. 



1 



— lim 



£„ = --iim j \{e)\w{e)\im[[w{e)]-^r 



w*{e) 

Eqns. (|32p and ([55)1 can be written as follows: 



N 



CnZ ^ ^ Cn 

n=0 n=N+l 







(34) 
(35) 

(36) 



N 
n=0 



Jm 



w{e) 



£nli-nlim[[Wi0)Y 

n=N+l 
= 



(37) 
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In Eq. ((36|) . we substitute for c„ from Eq. (|35p and solve 
for a so that, 



(z2) 



2^JV+l 



AT 



J - ^ c„z" 



?i=0 



(l^z^) 

(Z2) 



When no space-like constraints are used, Eq. ([33]), reduces 
to the A equation in [TTlfT^ . We note that [T3j obtains the 
A equation keeping only the charge radius fixed. 

We can further simplify the expression for a, given in 
Eq. dSHD, to get, 



. / A(.)[^(^j'(-)""--"?%^ 



r 



W*{e) l-re- 



(38) 



(z2)W+i 



^ AT+l 

J - y c„z" -t- / d0'x{0')p{e') 



Next, we simplify the equation for time-like constraints 
(Eq. ([37]) ). using c„, given in Eq. ([35]) . so that we have. 



(46) 

The simple condition = J^o'i^^n)'^ < 1 [H] now gets 
modified as follows 



N 



c„Im 



W^(6') 



ay z hm Im 

n=Ar+l 

- 1 limlm / A(0')^^^T^ y r"e™(^-^')d0' 
Trr^i Jr ^ 'w*{d')W{e) 

= 0. (39) 
Using the following: 

W{e) = |VF(0)|e*'^(^) 

where, 



n=0 



- E (2")'<1- (47) 

n=N+l 



We can once again substitute for c„ and use Plemelj- 
Privalov relations (|43|) to simplify the algebra so that. 



(40) 



= E(cn)' + ^ E c„ /" d(?A(0) sin H - ^{0)] 

n=0 n=0 "^-^ 

\ n=0 / 



(41) 



In our case, we will set N = 2> following TT] . The region 
r is defined as follows: 



and 0(0) is the phase of the outer function W7r(C = e*^), we 
substitute for W{6) in Eq. ([351) ■ that can be subsequently 
simplified to the following: 



r = 



{0, 0i„} 
{27r-6'i„,27r} 



(49) 



2^,"gi("e-<7!.(e)) 



> c„ sm [n6l - <P{e)] + a lim ) Im ( 



(42) 



Using Plemelj-Privalov relations [TT] , 



1 ^.^ /■^'^ ^^(6*') 



1 

^(^^) + -V - 

71- Jo 1 



F(0') 



and the expression for a, given in Eq. psp . we can simplify 
Eq. ([31]) to get the following equation for X{9). 



N 



= - x{e) + E 



n=0 



sin(n0-<Z>(^^))--^/3(0)z 



Tri^ ^ ^2 sin[^] 
+ i 1^ d0'A(0')(l - ^')/3(^)/3(e') + J^^PiS) (44) 
where 

_ sin [{N + l)g - g'(g)] - z sin [iVg - g>(g)] 
^ l + z2_2zcos(0) ^ ^ 



From Eq. ([?7)) . it is clear that information on the phase 
of the form factor along different parts of the cut can be 
used, i.e.. Si can be piece- wise continuous. 

We solve for X{9) using eqns. ([33]) ■ obtain the corre- 
sponding value of a from Eq. (|46]) . Using c„'s, already 
constrained by normalization of -FV and pion charge ra- 
dius, fiQ is evaluated (Eq. (|48p ). Only those coefficients, 
c,i, which satisfy //g < 1 are retained. 

We have carried out a thorough check of our results by 
comparing with similar results in the literature. In partic- 
ular, our results for the case of = 2 can be readily 
checked against results in [H] , that were obtained for the 
case of TT X scattering, where the Callan-Treiman point 
was used. Although time-like, the Callan-Treiman point 
lies within the region of analyticity, i.e., \z\ < 1, and as a 
result, can be introduced through linear constraint equa- 
tions. 



3 Data 

In this section, we present the data used to constrain c 
and d. We use space-like data from recent experiments [U 
[9] . These are summarized in [2] . We present the data here 
once again for the sake of completeness. 

The time-like phase we will use for part of this study 
is defined as 



Sl{t) = arc tan 



■ni?p — t 



(50) 
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Table 1. Space-like data from Tadevosyan et al. [9] and Amendolia et al. [8], used in this work. 





i(-Q2) [GeV^] 




zit) 


h{z) X 10"^ 


Tadevosyan 


-0.600 


0.433 


± 


0.017 


-0.494 


2.915 


Amendolia 


-0.131 


0.807 


± 


0.015 


-0.242 


4.454 




Fig. 2. Percent discrepancy between the phases obtained from 
the analytical model used in [TT] {Sp) and those from the Roy 
equation using two different sets for the scattering lengths [161 
[17]. 



^ 



— Cousu-alaui froir. VJS') = ' ant 

— timelike phase 




Fig. 3. Results for constraining c and d using the normaliza- 
tion -F,r(0) = 1 and the Pion charge radius, r-^ - large ellipse. We 
also see (small ellipse), the bounds obtained when the phase of 
the time-like phase alone is used in addition, using as input. 



and 

where nip = 770 McV is the mass of the p meson, Fp — 
150 MeV is the width of the p resonance and 771,^ = 139 McV 
is the mass of the Pion. This parametrization first pro- 
posed in ref. (TS] and adopted in ref. |TT] has been used 
here for ease of comparison of our results with those in the 
literature. It may be noted that at low energies, Eq. ((50|) . 
agrees well with the one-loop chiral perturbation theory 
expression for the two-pion elastic scattering phase shifts 
and also with experiments for t > 0.5 GeV^, as noted 
in [11]. Therefore, we assume that the phase of the pion 
form factor coincides with Eq. ([50]) for t^^ < t < tin, where 
iin = 0.8 GeV^. Our conformal map gives the following 
relation between t and 9 for z = C = e*^ 

t ^ U + UtW{e/2) (52) 

from which we can calculate ^in for ijn- 

In order to test the sensitivity to the parametriza- 
tion of the phase, we adopt the phase shift obtained from 
Eq.(Dl) of the accurate Roy equation study of Ref. [16] . 
for which the inputs are taken to be (aQ,aQ). Interest- 
ingly with the choice (a[j,a§) = (0.225, -0.0371) [16] the 
phase shifts do not differ from this analytical model used 
in [11| by more than 5.5% in the entire range, where as 
with the more recent favored value of (aQ,a§) = (0.220, 
—0.0444) [17], the difference can be as much as 10% as 
shown in Fig. 0. Despite this, the influence on our final 
results is not appreciable as we will show. 



The Pion form factors can be directly determined by 
the Kiihn-Santamaria and Gounaris-Sakurai parametriza- 
tion [IH] . We can, alternatively, use these two parametriza- 
tion for the form factor, obtain the phase and check the 
sensitivity of the results obtained using Eq. ([50]) . The cen- 
tral value of the parameters given in |18] are used for the 
fits we consider in this work. As we will see, the general 
lack of sensitivity to the parametrization of the form fac- 
tor, implies that no significant improvement will result 
from more recent or more precise inputs. 

In the next section we present our results for the above 
choice of parameters and data and also test the sensitivity 
of the bounds to variations in the data, choice of phase etc. 



4 Results 

In this section, we present our results for the bounds on 
the expansion coefficients and also check the sensitivity of 
the bounds to the errors in the input information. 

Fig. ([3]) shows the bounds obtained using just the nor- 
malization of and the pion charge radius (large 
ellipse). The smaller ellipse, shows the improvements on 
the bounds when the phase of time-like data, as given in 
Eq. ([50]) . introduced through the constraint Eq. ([28]) . is 
used. The constraints obtained on the Taylor coefRcients 
can be further improved by including the space-like data 
along with the phase of time-like data, discussed in sec- 
tion as seen in Fig. ([!]). The open ellipse is the allowed 
region in the c-d plane when one space-like datum, from 
Tadevosyan et al (refer Table [I}, is used and the filled 
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= time-like phase 
500 - — dme-lilce phase + 1 space-hke 
® 1 space-like 



> 

o 




c [GeV 



Fig. 4. Results for constraining c and d using pure time- 
like (squares), pure space- like (circles) and combined time- like 
phase and one space-like datum (solid-line) from Tadevosyan 
et al (see Table [1]). The insert shows the area of intersection 
where the combined bounds lie. 




c [GeV 1 



Fig. 5. Results for constraining c and d using phase of time- 
like data along with one space- like data (value closest tot = 0), 
for the sets [8]|9], as in table 
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Fig. 6. Variations in the bound as the space-like data from 
table (m are varied within the error bounds. 
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Fig. 7. 

changed 
Eq. dMl) 



c [GeV" ] c [GeV" ] 

Variation in the bounds as the phase information is 
5p is the parametrization for the 5l{t) defined in 



one represents that obtained when pure time-hke phase 
is used. Now when both phase of time-like and a single 
space-like datum are combined, the allowed region is an 
ellipse in the region of intersection of the respective el- 
lipses (refer inset in Fig. ([4])). As is evident, the overlap 
region of the pure space-like and pure phase of time-hke 
is significantly larger than the true region where they are 
taken simultaneously. 

We now focus on the bounds obtained when time-like 
phase and space-like datum are simultaneously used and 
address the sensitivity of these bounds to the variations 
in the input. Fig. ([5|) shows the variation in the bounds 
obtained as the space-like datum moves away from t = 0. 
The data point h{z) from Amendoha, given in table ([!]), is 
closer to t = compared to the corresponding data point 
from Tadevosyan (table ([1]). The figure clearly shows that 
space-like data from low t value constrain the coefficients 
better compared to data from higher t value. We note 
that the same trend was observed for pure space-hke con- 
straints [2]. 



Fig. ^ shows the variations in the bounds when the 
space-like datum in table H]) is varied away from its cen- 
tral value. Here, "max" refers to the central value -I- maxi- 
mum error and "min" is the central value — maximum er- 
ror. We see that the data is indeed sensitive to the errors 
in the space- like datum. However, varying the time-like 
phase between the definition in Eq. (jSO]) and direct eval- 
uation from Kiihn-Santamaria (KS) or Gounaris-Sakurai 
(GS) fits for the form factor docs not change the bounds 
on c and d, as seen in Fig. ([7]). This can be attributed to 
the fact that the phases produced by GS and KS fits and 
Eq. (|5D|) agree with each other for the range of t considered 
here. GS and KS parametrization involve two additional 
resonances compared to Eq. ([50|) that lie outside the range 
of t used here. 

We have also carried out an analysis with accurate Roy 
equation fits to the phase shifts. From Fig. ([8]), we observe 
that the effect on the bounds of the coefficients c and d 
is rather weak as noted for other fits like KS or GS. The 
value of tin, here fixed at tin = 0.8 GcV^, defines the inter- 
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Table 2. Bounds from phase of time- like data for different observables 





n' 11 




d[GeV-^] 
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Fig. 8. Comparing the constraints obtained when the fits given 
in Eq. ([50} are compared to the more accurate ones from the 
Roy equation |16II17] . Note that the dependence is negligible. 
The space-like datum comes from the set of Tadevosyan et al. 
(Table [1]). 
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Fig. 9. Varying the integration region for the phases ob- 
tained from the Roy equation |16lll7j . The left panel uses 
tin = 0.8 GeV'^, while the right panel uses tin = 0.6 GeV^. 
Note the dependence on tin, as well as on the two different pa- 
rameter choices for the phase shifts from the Roy equation, is 
weak. 



val r, as given in Eq. (|^5)) . Strictly speaking the solution 
(Eq.(DT) of Ref [IS]) is valid up to t^ax of 0.64 GeV^ but 
we assume that it is valid until tniax = 0.8 GeV^, as there 
is good agreement with the model phase Sp (see Fig. [2]). 
Furthermore, as shown in Fig. dH), lowering tin to 0.6 GeV^ 
does not lead to a perceptible change in the allowed region 
for (c, d) . Similarly, varying tin for the analytical model dp 




Fig. 10. Variation in the bounds as the interval F is varied. 
The phase is obtained from the Sp model given in Eq. 1501 By 
changing tin, more or less of the phase information is included. 



for the data sets from Amendolia and Tadevosyan has a 
very weak effect on the bounds, as in Fig. HI 



5 Discussion and Summary 

In this paper, we study the improvements on the bounds 
of the low-energy Taylor expansion coefficients of the Pion 
EM form factor, when both the phase of the time- like data 
and one space-like datum is used. We use the method of 
Lagrange multipliers to include the constraints. Bounds 
obtained using just the phase of time-like data with ei- 
ther or 77' jTl] as input is shown in table ^ The re- 
sults, when time-like phase and a single space-like datum 
are simultaneously used, are encouraging. In this work we 
have carefully considered the effect of a variety of ways in 
which the time-like phase is accounted for, by considering 
a p dominant model, phases from GS and KS parametriza- 
tion of experimental data and accurate Roy equation fits 
to scattering phase and have demonstrated that the results 
are stable. Our best estimates are obtained when space- 
like data from the data set of Amendolia at al. [5], as given 
in table (H]), is used. In this case, the coefficient c lies in 
the range: [2.3 GeV-^ 5.4 GeV""] and d in [-14GeV-^ 
56GeV~^]. Furthermore, if the modulus of the form fac- 
tor is included along with the phase [TT] and U' as input, 
the c and d have the following ranges: [0.5 GeV^, 7.5 
GeV"^] and [-1 GeV~^ 22 GeV"^]. It is worth noting 
here, that the range of c isolated in the present study, is 
the most stringent and is also in agreement with chiral 
perturbation theory, where c has been determined to 2- 
loop accuracy: 4.49 GeV~^ [12] (see also ref. [lO])- As an 
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interesting exercise, fixing the value of c to that obtained 
from chiral perturbation theory, we read-off the range for 
d. Both time- like phase and time- like phase plus one space- 
like constraint gives d in the range ~ 20 — 30GcV^®, 
when ttfj, is used as input. Similar range is obtained for 
the bounds on d when U' and phase of time-like data are 
used as inputs. When the modulus of the form factor in 
the time- like region is also taken into account, along with 
the time-like phase information, the range of d is roughly 
around ~ 10 — 15 GeV~^. The slight mismatch between 
these ranges could warrant a separate study taking into 
account the uncertainties in parametrization of the phase 
and modulus and corrections to U' and is clearly beyond 
the scope of this work. However it is remarkable that such 
a different variety of inputs and theory leads to a rather 
coherent picture for the values of c and d. Finally it may 
be noted that the value of the bound on d is an order of 
magnitude greater in this case, when pure space-like data 
is considered [5] 

The GS fit is well known for its good analytic proper- 
ties. Therefore, the value of c and d can be obtained by ex- 
panding out the parametrization for the form factor. This 
gives an estimate for c,d as 3.37 GeV"'' and 10.2 GeV~^ 
respectively and is well accommodated by our best con- 
straints obtained for the Amendolia data set for the space- 
like part. The value of c so determined is consistent with 
another determination available in the literature ref. [21j 
of 3.2 ± 0.5 ± 0.9 GeV"''. The more accurate fits for the 
phase shifts from the Roy equation [TBlflT] yield bounds 
which arc very close to the ones obtained from the ana- 
lytical model 6 p. Fixing the upper limit of the integration 
region F to be tin = 0.8 GeV^ and the space-like datum 
to that from Tadevosyan et al. (table [T]), the bounds on 
c from the model Sp for the phase shift is: [—0.36, 11.4], 
which is identical to that obtained from the Roy equation 
fits using parameters from [T6jll7j. The bounds obtained 
for d varies slightly for different phase shift parameteriza- 
tions, i.e., [—51.12, 142.51] for the model phase shifts and 
[—55.73, 142.44] for the Roy equation fits using parame- 
ters from [T7|; while the range obtained for d using the 
parameters from |16| is identical to that obtained from 
the model phase shifts. From this we can conclude that 
the Roy solutions and the model Sp essentially agree in 
the region of integration, especially in the resonance re- 
gion. 

At this point, it becomes important to include the 
modulus and phase of time-like data, along with space- 
like data, so that we include all the information available 
and constrain the pion form factor using arguments of 
analyticity. Hence, it would be interesting to explore the 
possibility of working along the lines of [11], where the 
modulus and phase of time-like data are used to obtain 
bounds, and combine it with the technique already ex- 
plored in [2j. 

We note that the bounds are sensitive to the errors in 
the data used. Therefore, as an intermediate step, it would 
be worth-while to obtain a theoretical fit to the data (x^ 
fit) and use the fit as input. This would eliminate the in- 
fluence of experimental errors. In order to completely un- 



derstand the sensitivity issues, it is important, as already 
noted in pjl22j. to do a combined error analysis of time- 
like and space-like data, based on the work of Raina and 
Singh [22] , that also gives the possibility of including the 
modulus via the technique of Lagrange multipliers. This 
theory needs to be developed. 
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